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Quantum Dimer Models (QDM) arise as low energy effective models for frustrated magnets. 
Some of these models have proven successful in generating a scenario for exotic spin liquid phases 
with deconfined spinous. Doping, i.e. the introduction of mobile holes, has been considered within 
the QDM framework and partially studied. A fundamental issue is the possible existence of a 
superconducting phase in such systems and its properties. For this purpose, the question of the 
statistics of the mobile holes (or "holons") shall be addressed first. Such issues are studied in details 
in this paper for generic doped QDM defined on the most common two-dimensional lattices (square, 
triangular, honeycomb, kagome,...) and involving general resonant loops. We prove a general 
"statistical transmutation" symmetry of such doped QDM by using composite operators of dimers 
and holes. This exact transformation enables to define duality equivalence classes (or families) of 
doped QDM, and provides the analytic framework to analyze dynamical statistical transmutations. 
We discuss various possible superconducting phases of the system. In particular, the possibility of 
an exotic superconducting phase originating from the condensation of (bosonic) charge-e holons is 
examined. A numerical evidence of such a superconducting phase is presented in the case of the 
triangular lattice, by introducing a novel gauge- invariant holon Green's function. We also make the 
connection with a Bose-Hubbard model on the kagome lattice which gives rise, as an effective model 
in the limit of strong interactions, to a doped QDM on the triangular lattice. 



I. INTRODUCTION 

In 1987 Anderson^ suggested the strange behavior of 
cuprate materials between the superconducting dome 
and the magnetically ordered insulating phase could be 
described by a Resonating Valence Bond (RVB) state 
in which preexisting magnetic singlet pairs of the insu- 
lating state become charged superconducting pairs when 
the insulator is doped. Just one year later appeared the 
first effective model in which the magnetic degrees of free- 
dom are disregarded in favor of the more pertinent singlet 
degrees of freedom^. This is nothing but the Quantum 
Dimer Model (QDM). However, it was soon realized^ that 
the groundstate of the undoped S* = 1/2 Heisenberg an- 
tiferromagnet on a square lattice exhibits a long-range 
Neel order in contrast to the initial expectation based on 
the RVB picture. Furthermore, the QDM on a square 
lattice was also found to have only gapped crystalline 
phases but no evidence of an RVB spin liquid phase in 
a finite region of the phase diagram.^ It is still possible 
to argue that, even though the undoped antiferromag- 
net has the Neel order, the RVB picture gives a better 
theoretical starting point once the system is doped with 
holes. However, it would be natural to ask if the RVB 
spin liquid phase can be realized in undoped magnetic 
system with only short-range interaction. 

One may expect that magnetic frustration would favor 
RVB state over the Neel phase. Thus, over the time, the 
main interest in QDM was shifted from the original moti- 
vation of the application to high Tc superconductivity, to 
the effects of frustration. However, clear confirmation of 
a RVB phase remained elusive for a rather long time. A 
breakthrough in the study of QDM was due to Moessner 
and Sondhi^ who showed that a simple QDM defined in 



the triangular lattice exhibit a disordered phase which, 
recalling that these dimer models are supposed to be ef- 
fective models for frustrated magnets, can be considered 
as an explicit example of the RVB spin liquid phase. It 
was also recognized that, the RVB spin liquid phase is a 
topologically ordered phase with a nontrivial topological 
degeneracy of the groundstates.^ In fact, the RVB spin 
liquid phase is essentially identical to the Z2 topologi- 
cal phase which was introduced in a completely different 
context of quantum information processing.^ The QDM 
is generally not exactly cqivalcnt to an antiferromagnetic 
Hamiltonian defined in terms of quantum spins. How- 
ever, the projection from a magnetic system to a QDM 
was performed successfully in Heisenberg antiferromag- 
nets on frustrated lattice, such as the square lattice with 
strong enough second and/or third neighbors couplings^'^ 
or the kagome lattice^°. These suggest that the QDM 
may well represent phases without magnetic order in an- 
tiferromagnets. In fact, very recently, frustrated Heisen- 
berg antiferromagnets on kagome and other lattices are 
reported to be in the RVB spin liquid phase (Z2 topolog- 
ical phase) by several authors. ^^"^"^ 

Now that the existence of the RVB spin liquid phase 
appears to be confirmed in QDMs as well as in anti- 
ferromagnets, the issue of superconductivity in doped 
spin liquids becomes a more pressing question. This iss- 
sue started in fact to be investigated shortly after the 
apearence of QDM^^ Doping of an RVB spin liquid is ex- 
pected to induce a novel type of elementary excitations 
called holon. A holon, carrying electric charge e but no 
spin, appears as a result of fractionalization, namely de- 
confinement of fractionalized excitations. Indeed, topo- 
logical degeneracy of the undoped RVB spin liquid is 
known to be intimately connected to the fractionaliza- 
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tion phenomenon^^ . 

Superconductivity may be realized if the holons con- 
dense. At least naively, one may expect that the resulting 
superconductor is an exotic one due to condensation of 
charge-e holons, instead of usual chargc-2e Cooper pairs. 
A fundamental issue in this problem is the statistics of 
the holon. For the holons to condense without form- 
ing pairs, they must be bosons. However, it should be 
noted that transmutation of the statistics^^'^^ is possible. 
Namely, the statistics of holons as elementary excitations 
appearing in the low-energy limit can be different from 
the statistics assigned to holons in the microscopic model. 

In this paper we address the issue of the statistics 
of holes and its interplay with possible superconducting 
phases in doped QDMs. In a recent work^'^ it was shown 
that a QDM with fcrmionic (at microscopic level) holes 
is equivalent to another QDM with bosonic holes. Be- 
cause of the equivalence, the statistics of the holon as 
a physical elementary excitation must be the same for 
either representation. This proves the existence of a dy- 
namical statistical transmutation in the system. In this 
paper wc study in more details the statistical transmu- 
tation in QDM and give a simple and efficient method 
to obtain the relation between the QDMs with fcrmionic 
and bosonic representation of the holes. 

In section II wc introduce a second quantization nota- 
tion for QDM Hamiltonians and show the gauge symme- 
try associated with them. In section III we present the 
composite particle representation of QDM Hamiltonians 
which is the key ingredient to show the exact equivalence 
between a QDM with bosonic and another QDM with 
fcrmionic holes. This equivalence is shown for a generic 
flipping term defined in any kind of lattices. The result, 
which relies in an orientation prescription of the bonds in 
the lattice considered, is totally generic and can then be 
applied to any QDM defined in the most common lattices. 
It is important to stress that the method used here differs 
considerably with the one used in Ref. 18 where a two- 
dimensional version of the Jordan- Wigner (JW) trans- 
formation was used. As we explain in the conclusion, 
it is important to point out the numerous advantages of 
this new representation compared to the JW transfor- 
mation. In section IV we argue how the modification of 
the orientation prescription can be interpreted as a sim- 
ple gauge transformation in the QDM Hamiltonian. We 
then apply the general result of the statistical transmuta- 
tion obtained in section HI to generic QDM Hamiltonians 
defined on the square, triangular, hexagonal and kagome 
lattices. Section V is devoted to numerical investigation 
of four inequivalent QDM defined on the triangular lat- 
tice. In particular, we identify an exotic superconductor 
phase due to condensation of holons with charge e, mea- 
suring the gauge-invariant Green's function of a single 
holon. In section VI we discuss an explicit realization of 
one of the QDM discussed in section V. It is obtained as 
a low energy strong interaction limit of a Bose-Hubbard 
model on the kagome lattice. The number of bosons is 
directly related to the doping, or number of holes, in the 



resulting QDM on the triangular lattice. Section VII is 
devoted to the discussion of our results. We also include 
as an appendix the derivation of the statistical transmu- 
tation for a generic QDM on the kagome lattice using 
the Jordan- Wigner transformation. Of course the result 
is consistent with the one obtained with the composite 
particle representation obtained in section HI, but allows 
a better understanding of the connection between these 
two different methods. 



II. THE HAMILTONIAN AND ITS GAUGE 
SYMMETRIES. 

We start with a doped quantum dimer model on a two- 
dimensional lattice. To fix the ideas, wc work here with 
the Hamiltonian defined on the square lattice but all the 
arguments remain valid for any two dimensional lattice. 
We write the Hamiltonian as 



H = Hj + Hv + fl", 



(1) 



with 
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where the sums are over all the smallest resonant pla- 
qucttes on the lattice (for the square lattice these are 
the squares). In a second quantized formalism we as- 
sume that dimer configurations are created by spatially 
symmetric dimer operators bl ■ and holes are created by 

bosonic operators a\. Then, we can re- write the Hamil- 
tonian as: 



Hv 
Ht 



jT.[^iAM^i^^+^-^] (2) 

□ 

□ 

^E{4j^J'fc4a^ + H.c}. (4) 



In the last equation, the indices correspond to the label- 
ing of the sites of a square plaquette as in Figure 2. In 
our previous conventions, dimer configurations are rep- 
resented by spatially symmetric operators h\ j satisfying: 



h 



(5) 



The boson operator aj creates a hole in the site i and 
satisfies 



(6) 



ala] 



= 0. 



The operators a and b commutes one with each other 



[a,;, bj 



alb] 



= 0. 



(7) 



We introduce in the model a constraint on the number of 
dimers and holes which warrants that at each site of the 
lattice there is either one and only one hole or one and 
only one dimer arriving to it: 



E 



bl,+,b,,,+z 



= 1. 



(8) 



Of course this constraint implies, among others, that the 
holes have to be considered as hard core bosons. It is im- 
portant to notice that the Hamiltonian has the following 
U{1) gauge symmetry: 

(9) 
(10) 



2j -> e*'^^' aj 



^ e 



where is an angle. This invariance can be exploited 
to prove the statistical transmutation symmetry in some 
two dimensional systems by means a Jordan- Wigner 
transformation on the holon operators^®. In the following 
we present an alternative description for the doped QDM 
using composite operators which allows us to understand 
in a different way the equivalence between a model with 
bosonic holes and one with fermionic holes. For doing 
this, we have first to make a choice of a given orientation 
prescription for the bonds in the lattice. 



III. THE "COMPOSITE" REPRESENTATION 
FOR THE QDM 

A. Composite particles 

In order to prove the equivalence between a QDM 
Hamiltonian with bosonic holes and another Hamiltonian 
with fermionic holes, we propose a different formulation 
for the QDM. This formulation is done in terms of com- 
posite particles by defining the operator 



=6,jaja]. (11) 

This operator destroys a dimer between sites i and j and 
creates two holes at the same sites. Let's call He the 
subspace of states that satisfies the constraint (8). For a 
given state G Tic we have that \ip) — Bi,j\ip) is also a 
vector in "He- 
One can easily notice that the operator Bi j is invariant 
with respect to the gauge transformation 

aj e'^'aj (12) 

bjM ^ e'^^^+^-^^k (13) 




FIG. 1. (color online) Schematic snapshot of a doped "dimer 
liquid". Each site is occupied by either a (single) dimer or a 
hole (empty site). 




b) 



FIG. 2. a) Indexes corresponding to each square plaquette in 
the Hamiltonians Hj and Hv- b) Indexes corresponding to 
each hopping process in Ht- 



This U{1) gauge symmetry was exploited in Ref.ll to 
represent doped QDM as theories in which a matter field 
is coupled to a gauge field. More importantly, one can 
check that within the subspace He, the set of operators 
{Bij} form a closed algebra similar to the one of {bij}. 
The Hamiltonian can entirely be written in terms of these 
Bi j operators making its gauge invariance a trivial state- 
ment. Its precise form is given by iJ = Hj + Hy + Ht 
with: 



HJ^-JJ2 {BhBliBj,kBi,, + H.c.} 



(14) 



Hv = V 



E {BhBli^^^^Bk, + Bl,BlB,^kBu}{l5) 



Ht = -tJ2{BlB,^k+ii.c} 



(16) 



It is evident that, within this formulation, the basic 
building blocks of the model are created by Bj^ which 
corresponds to composite particles of charge 2e. Namely, 
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the model is completely defined in terms of the con- 
stituent particle with charge 2e. This has several impor- 
tant consequences. In particular, the gauge invariance re- 
quires that the energy spectrum of the system on a torus, 
as a function of the magnetic flux $ through the "hole" 
of the torus, is invariant under $—!>$ + 7r/e. (For early 
discussions on the Tr/e-flux periodicity in the QDMs, see 
Refs. 19 and 20 and references therein.) This periodic- 
ity corresponds to the unit flux quantum for charge 2e 
objects. However, this does not necessarily mean that 
the physical elementary excitations of the system have 
minimum charge 2e.^^ The system can have a topological 
order which leads to fractionalization; elementary excita- 
tions can have fractions of the charge 2e of the constituent 
particle of the microscopic Hamiltonian. If the charge- 
e holons are deconfined as a result of fractionalization, 
they could condense to form an exotic superconductor. 

The apparent contradiction between the periodicity of 
the energy spectrum in n/e flux and the expected flux 
quantization in the unit of 27r/e in the condensate of 
charge-e holons is resolved by the existence of the topo- 
logical vortex excitation called vison. Insertion of the ir/e 
flux corresponds to trapping of a vison. Although the flux 
periodicity of the groundstate energy does not distinguish 
an exotic charge-e condensate from the usual supercon- 
ductor, an experimental detection scheme of the charge-e 
condensation, based on a "vortex memory effect" , was 
proposed. An actual experiment'^^ on the high- Tc su- 
perconductor did not flnd such a signature of charge-e 
condensation. Nevertheless, the exotic supercondutiv- 
ity due to condensation of charge-e objects is possible in 
principle, and is an interesting subject to pursue theo- 
retically and experimentally. Later in this paper, we will 
introduce a quantity which detects a charge-e condensa- 
tion, and study it numerically in several QDMs. 

B. Statistical transmutation symmetry 

The main advantage of the formulation in terms of 
composite operators presented above is that one can 
prove an equivalence between a Hamiltonian where holes 
are hard core bosons and another one where the holes 
are fermions. Let us consider a QDM hamiltonian with 
bosonic holes, where their creation and annihilation op- 
erators a| and aj satisfy bosonic commutation relations. 
Let us also consider another QDM hamiltonian with 
fermionic holes, created and annihilated by the set of 
operators and fj which now satisfy fermionic anti- 
commutation relations. We then build the respective 
composite operators: 

B^.l - KA4 (17) 

= Kj flf] (18) 

As for the operators B, all quantities in the rest of the 
paper with a tilde correspond to operators and coupling 
constants of the fermionic representation for the holes. 



Before proceeding, there is an important statement to 
make. Again, one can show that within the subspace 
^-Lc^ both set of operators {5,;.^} and {-Bi.j} form the 
same closed algebra of bosonic dimer operators. Another 
important point is that the deflnition of the composite 
operators in terms of fermions is more subtle because it 
is necessary to take a prescription for the orientation of 
the dimers (which determines the order of the fermions 
in the endpoints of each dimer). 

In the following we will call "even prescription" of a 
given plaquette an ordering prescription for the bonds 
such that all the bonds are oriented in a clockwise di- 
rection or an even number of bonds are oriented anti- 
clockwise. By contrast we call "odd prescription" the pre- 
scriptions obtained from the clockwise ordering by flip- 
ping an odd number of bonds. 

Notice that, since the resonance plaquettes containing 
N dimers have necessary 2A^ bonds, the anti-clockwise 
prescription(where all bonds are oriented anti-clockwise) 
is always an even prescription. 



a. Theorem: Given a resonant plaquette of arbitrary 
length with an even prescription for the bonds, then, for 
the kinetic term of the dimers in the plaquette, we have 
the equivalence: Hj(J, bosons) O Hj{— J, fermions). 
In other words, the resonance term of dimers in the 
plaquette is invariant under a simultaneous change of 
"statistics" of the holes in the system (i.e. bosonic into 
fermionic or vice-versa) and the sign of the dimer reso- 
nance loop amplitude J. 

Proof: 

Consider a resonance loop containing 2A^ sites {N 
dimers) numbered from 1 to 2A'^ in the clockwise direc- 
tion as in Fig. 3-a and Fig. 3-b. The kinetic Hamiltonian 
for dimers belonging this loop can be written in terms of 
bond operators bi_j as 

^j' ~ "^[^1,2^3,4^5,6 ■ ' ■ ^2Af-3,2Af-3^2Af-l,2Af] i^^) 
X [&2,3&4,5---&2A'--2,2Ar-l&2JV,l] +H.C (20) 

where, the index N indicates the number of dimers in the 
loop. 

Now, we add one dimer (2 sites) to the loop, obtaining 
a resonance loop with + 1 dimers. In this case, the 
Hamiltonian can be written as 

Hj — J L"l,2"3,4"5,6 "2Af-3,27V-2"2Af-l,2A'"Q,^J 
X [&2, 3^4,5 ■•■ &2iV-2,2Af-l&2Ar,Q^/3,l] + H.C 

Since all the operators are acting on different bonds they 
all commute and we can rearrange them in the following 
way 

n J — J ['^i^2"3,4 "2Ar-3,27V-3"27V-l,2ArJ 
X [^2,3^4,5 • ■ • ^2^-2,2^-1] ^2^,1 + H.C 
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FIG. 3. Elements used to prove the inductive step, a) and b) correspond to the two possible dimerization in a A'' dimers 
plaquette whereas c) and d) correspond to the plaquette with A*" + 1 dimers. 



r/A^+i - /[At /,t ...!,t lA 1 

ri j — J [U^ 2^3,4 "2W-3,2Ar-3"2Ar-l,27Vj 

X [62,3^4,5 ■ • ■ &2JV--2,2Ar-l]foa,/3^2N,a6/3,l + H.C 

or in a compact notation 

( n ^1-1.2. j ( n ^2j.2j + l j &2iV,l + H.C. 



N 



H 



N+l 



•^1114-1,2, 




Now, wc insert on the right of the Hamiltonian the string 
of operator S^q ~ T^i=i fill ^ where the index i corre- 
spond to the sites on the resonance loop. This operator 
acts as the identity operator on the sites belonging the 
loop because fif} = 1 in the absence of holes. We start 



with the Hamiltonian 

N 



IJ2N,1 



X A//--- W2V + H.C 

Now, we move the fermions to the left in order to form 
the composite operators Bi j = bij/^f^ corresponding to 



the dimer operators inside the brackets. Commutation of 
the fermions gives a global sign. 



N 



H 



N 




1 -y-^fJ\\[Bh-.. 

We can follow exactly the same procedure in the loop 

with + 1 dimers, the global sign resulting from the 

commutation of fermion operators to write the products 

in terms of composite particles is the same that in the N 

dimers case. We can write for the + 1 case 
H.C 



N 




X bl^pb2N,abp^iflfl^fafifpfl + H.C 

Now we can determine the change in the sign of J 
when a dimer is added in the loop. First wc commute 
the operators fl and flj^ in to form the operator 

B2N,i = b2N,ji.f2Nfi- This commutation gives another 
sign to complete the global phase in the Hamiltonian. 
Then wc can write for the A^ dimers case 

^ |^n^2j-i,2jj ^2,,2,+ijB2A.a + H.c 

In the Hamiltonian correspoirding to A/ + 1 dimers we 
have to define 3 composite operators. We have that 



= (-1) {^Ipfpfo) {b2N,afl^fi) {bp.lflfl 



( ^ 1 ) -^I , ;3 -^2 AT , Q -B^J , 1 



(21) 
(22) 
(23) 
(24) 



I 

Finally, the Hamiltonian corresponding to A'^ + 1 dimers reads 



N 



ffr'=(-i)«+v|^n^2,-M, 

^ BB2N,aBp^i + H.C 
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We have proved by induction that, if for a N dimers 
loop, the kinetic term acquires a given sign when the 
Hamiltonian is written in terms of fermionic composite 
particles, then the kinetic amplitude corresponding to 
A'^ + 1 dimers acquires the same sign. To complete this 
mathematical induction proof we need check that the 
statement holds for the lowest value of N. The small- 
est possible resonance loop is given by a loop with only 
2 dimers. It is easy to check that in this case 

Hj = Jb\,2bl,ih2,^hi^ + H.c 

= (-l)JBj 25^,4^2,354,1 + H.c 

Then we have proved that the kinetic term for a res- 
onance loop of arbitrary length oriented in a clockwise 
direction, the amplitude J in the Hamiltonian written 
using dimer operators bi^j changes to — J when we write 
the Hamiltonian in terms of fermionic composite opera- 
tors Bij . A trivial verification shows that the amplitude 
remains unchanged when we write the Hamiltonian in 
terms of bosonic operators Bi^j = bija\a}y 

The result above can be re- written in the more appeal- 
ing way: 

Hj{J,B)=Hj{~J,B) (25) 

The proof can easily be extended to the potential term 
Hy- In this case is easy to see that the bosonic and 
fermionic versions give the same sign in the amplitude 
V. 

Hv{V,B) = Hv(V,B) (26) 

The equivalence proved above is valid for any even pre- 
scription on the plaquette. Starting from the clockwise 
prescription where the results above has been proved, 
if we flip two bonds this induce the commutation of two 
fermionic operators and the sign remains unchanged. But 
if we flip an odd number of loops we must to commute an 
odd number of extra fermionic commutations in order to 
form the composite operators. These permutations gives 
an extra sign in the Hamiltonian. Then is easy to prove 
the following corollary: 

Corollary 1: 

Given a resonant plaquette of arbitrary length with an odd 
prescription for the bonds, then, for the kinetic energy 
of the dimers in the plaquette, we have the equivalence: 
Hj{J, bosons) O Hj{J, fermions). 

The equivalence in the potential term do not change if 
we take an odd or even prescription. Using this property 
of the potential term and the Corollary 1 we can derive 
the following corollary. 

Corollary 2: 

Note that we have actually proved that, if in a given lat- 
tice we can take an even prescription for all the plaquettes 



involved in the Hamiltonian. Then the equivalence 

H{J, y, bosons) = H{—J, V, fermions) (27) 

is valid for the Hamiltonian in the complete lattice, 
whereas if we can take an odd prescription for all the 
plaquettes in the Hamiltonian. We have the equivalence 

H{J, V, bosons) = H{J, V, fermions) (28) 

In order to complete the panorama for the doped QDM 
we study the fermionic and bosonic representation of the 
Hamiltonian Ht corresponding to the hopping of holes. 

Consider 3 nearest-neighbors sites of the lattice as in 
Fig. 2-b. In terms of bosonic holes and dimer operators 
we can write a general hoping term as 

hS,k = blb.^kaia, (29) 

if there is no hole in the intermediate site j we can add 
on the right the identity as aja| = 1. We then obtain 

^tik = -tbljbj^kaia.aja] (30) 
= {bj,ka]al^ (31) 

^-tB^B.^k (32) 

In the bosonic case we don't need to worry about the 
prescription in the lattice but it is important when we 
study the fermionic description. In this case we take the 
prescription i — > j — > fc. Starting from the Hamiltonian 

= blb,,,flf, (33) 

we insert on the right the operators fjfj = 1 

hS,k--tbl.,hkflhfjf^ (34) 

= t {by,.U) {hkfUl) (35) 

= -tBlB,., (36) 

Using the prescription i j — >■ k, the amplitude in 
the hopping term for the holes is the same if we use the 
fermionic or bosonic versions of the composite operators. 
Flipping two arrows we have the prescription fc — > j — i. 
It is a simple matter to see that with this prescription 
the hopping amplitude is also the same for the two cases. 

Then, the hopping of the holes written in terms of 
bosonic and fermionic composite operators have the same 
amplitude t, provided that we use one of the two prescrip- 
tions satisfying that the intermediate site has one incom- 
ing and one out-coming arrow. If we take this i ^ j ^ k 
prescription in all the sites of the lattice the arrows follow 
a sort of Kirchhoff's first rule: see Fig. 4- (a), (b), (d). 
We will call this kind of prescriptions as "zero-current" 
prescriptions. Of course it is only possible to satisfy this 
prescription in all the sites if the coordination number of 
the lattice is even. An example where this is not possible 
is the Honeycomb lattice (with z = 3). In this lattice 
it turns out that it is not possible to take a prescrip- 
tion with the same number of incoming and out-coming 
arrows in each site. See Fig. 4-(c). 



FIG. 4. (Color online) Bond prescriptions on the square (a), triangular (b), honeycomb (c) and kagome (d) lattices. Light-blue 
plaquettes have an even prescription while the green plaquette has an odd prescription. 



IV. QDM CLASSIFICATION FOR DIFFERENT 
LATTICES 

A. On the choice of the bond orientation 
prescription 

As we saw in the last seetion, in order to prove the 
equivalenee between Hamiltonians built with bosonic and 
fermionic operators, one needs a bond orientation pre- 
scription for the fermionie case. Of course, this prescrip- 
tion is totally arbitrary and before proceeding it is impor- 
tant to clarify the issue of a different choice of prescrip- 
tion. Let us imagine a generic lattice for which we have 
chosen two different prescriptions, A and B. To clarify 
the ideas, imagine that the orientation of all the bonds 
in prescription B are the same that in prescription A, 
except for one single bond, which is connecting points i 
and j. Then, starting from a bosonic Hamiltonian, by 
doing the transmutation, we end up with two different 
Hamiltonians Ha and Hb which have the same signs for 
all the flipping and hopping terms except for the ones 
that contain the bond ij. Let us illustrate this with the 
following example; consider the square lattice in which 
prescription A is the one given in Fig. 4. Then, imagine a 
prescription B where only the arrow between sites i and 



j is reversed, as shown in Fig. 5. Starting from the same 
bosonic Hamiltonian, after the statistical transmutation, 
we get the Hamiltonians Ha and Hb- What is the differ- 
ence between Ha and /^s? They have the same signs for 
all the flipping terms, except the the ones of plaquettes 
a and /3, which are the only two containing this reversed 
bond. Also, all hoping terms are the same except those 
containing the link ij. 

Although one could naively think that these two re- 
sulting Hamiltonians are not equivalent, in fact they are, 
as can be easily seen by performing the following gauge 
transformation: 

&n,m -> -^«,m itn^i, m = j or n = j, rn = i 
bi,m bi n, else. 

Most generally, it is easy to convince oneself that different 
choices of prescriptions give rise to apparently different 
Hamiltonian which in fact are equivalent under a certain 
gauge transformation. We are now going to consider each 
lattice in detail and justify for each of them the choice of 
prescription we have made. 
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A B 

FIG. 5. The change of prescription corresponding to revers- 
ing the orientation of one single bond {ij in the figure) corre- 
sponds to a gauge transformation where only configurations 
containing a dimer in the ij bond have their sign changed. 
This in turn has the effect of reversing the sign of the flipping 
and hoping terms containing the bond ij, as for example the 
flipping of plaquettes a and /?. 

B. Square lattice 

For the square lattice, we consider the prescription 
given in Fig. 4-(a). Using this prescription, the hopping 
amplitude (t) remains equal for the bosonic and fermionic 
representation B and B. On the other hand, the ki- 
netic amplitude corresponding to dimers (Ja) changes its 
sign if an even prescription is induced in the plaquette 
of length a. In Fig. 4 the plaquettes of lowest order are 
shown. Light-blue areas correspond to even prescriptions 
induced in the plaquettes while green areas correspond 
to odd prescriptions. The relative sign between the cou- 
plings in the fermionic and bosonic representations cor- 
responding to the 8 smallest plaquettes are presented in 
table I 



N 


Loop 


j/i 
j/t 


2 


u 


-1 


3 


LU 


1 


4 


1 1 1 1 


-1 






-1 






-1 


5 


1 


1 


1 






1 


1 








1 


1 


1 



TABLE I. Values of Jc,/ Ja corresponding to the lowest orders 
of the resonant plaquettes on the square lattice. 



have the same sign. Again, the change in the sign when 
we change from a bosonic representation of the holes to a 
fermionic one is determined by the parity of each flipping 
term. In Table II we show the results for flipping loops 
containing up to three dimers. 



iV 


Loop 


j/t 
j/t 




\A 


-1 


2 




-1 






-1 






1 


3 


f 


1 






-1 



TABLE II. Values of Ja/Ja corresponding to the lowest orders 
of the resonant plaquettes on the triangular lattice. 



D. Honeycomb lattice 

The case of the Honeycomb lattice is more subtle. The 
coordination number in this lattice is z = 3 and it is not 
possible to take a "zero-current" prescription. Therefore, 
it is not possible to find a prescription in which all the 
hopping terms would remain the same after the trans- 
mutation. We then use the prescription of Fig. 4-(c) in 
which all the hopping amplitudes for the holes change 
the sign when we change to the fermionic representation 
of the operators (i„ = —tn)- 

The relative signs between the ratios J/t and J/t are 
presented in table (III) for plaquettes of 3, 5 and 6 dimers. 



N 


Loop 


J/i 
J/t 


3 




-1 


5 




-1 


6 




1 



TABLE III. Values of Ja/Ja corresponding to the lowest or- 
ders of the resonant plaquettes on the honeycomb lattice. 



C. Triangular lattice 

For the triangular lattice, as in the square case, the 
coordination number is even and we can take a "zero- 
current" prescription as shown in Fig. 4- (b). Then 
the hopping amplitudes for bosonic and fermionic holons 



E. Kagome lattice 

For the kagome lattice, we have chosen the prescription 
depicted in figure 4-(d). All the possible allowed fiipping 
terms up to 12 bonds arc depicted in the table (VI) of the 
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appendix. Is it interesting to note that these loops are 
all without exception even. This feature is not specific to 
loops of short lengths and one can convince oneself that 
all allowed flipping loops of arbitrary length are even. 
More over, our prescription choice is such that the hop- 
ing terms within one triangle remain invariant under the 
statistical transmutation. From this, one can conclude 
that a Hamiltonian with flipping terms {J;} and bosonic 
holes is equivalent to a Hamiltonian with fermionic holes 
and with the signs of all flipping terms reversed {—Ji}- 

One last remark one can make about the kagome lat- 
tice relies on its intrinsic flexibility. Take any triangle of 
it and change the orientation of the three bonds belong- 
ing to that triangle only. It is easy to see that with 
the new prescription all the hopping terms, including 
those belonging to the chosen triangle, do not change 
signs. However, the flipping terms containing one (and 
only one) bond belonging to that triangle will change 
their signs, i.e., these flipping terms in the transmutated 
Hamiltonian have the same sign as in the bosonic model. 
What this means is that, in contrast to the other lattice 
studied here, it is possible on the kagome lattice to build 
gauge transformations which leave invariant the sign of 
all hoping amplitude while changing the sign of "some" 
flipping terms (even locally). 



F. Example of application of the "statistical 
transmutation" symmetry 

To illustrate the power and extend of our results, we 
concentrate on a couple of concrete examples taken on the 
square and triangular lattices respectively. Let us con- 
sider the QDM defined on the square lattice with only 
two and three dimer flipping terms. These terms corre- 
spond to the first and second row of table (I). Its sibling 
model can be defined on the triangular lattice by just 
considering the terms with N = 2 and only the second of 
the terms with = 3 in table (II) . In principle we would 
have 16 inequivalent Hamiltonian on each case. However, 
our statistical transmutation result tells us that the num- 
ber of inequivalent Hamiltonians is smaller. Indeed, we 
find only 8 inequivalent Hamiltonians for the case of the 
triangular lattice, which we dubbed lo-, 11^, III<t and TVa 
where a — ± corresponds to the sign of the hoping term. 
From these 8 classes only 4 classes are inequivalent for the 
case of the square lattice. The smaller number of equiva- 
lence classes in the later case is due to to the equivalence 
i ^ — < which is valid for the square lattice but not for 
the triangular^®. The result is summarized in table (IV). 



G. Transformation of assisted terms 

As QDM have to be taken as low energy effective mod- 
els of frustrated anti-ferromagnets, it is important to see 
if other kind of term, apart from those already mentioned 
here, arise in the effective QDM Hamiltonians. Examples 



Statistics 


J4 




t 


Family 


Rcfs (□) 


Refs (A) 


bosons 


+ 


+ 


+ 


1+ 


26, 27, 28, 29 


18, 27, 29 


bosons 


+ 


+ 


- 


I 


26, 28 


18 


bosons 


+ 


- 


+ 


11+ 






bosons 


+ 


- 


- 


II_ 






bosons 


- 


+ 


+ 


III+ 


26 


18 


bosons 


- 


+ 


- 


III- 


26 


18 


bosons 


- 


- 


+ 


IV+ 






bosons 


- 


- 


- 


IV- 






fcrmions 


+ 


+ 


+ 


III+ 


26 


18 


fcrmions 


+ 


+ 


- 


III- 


26 


18 


fcrmions 


+ 


- 


+ 


IV+ 






fcrmions 


+ 


- 


- 


IV- 






fcrmions 




+ 


+ 


1+ 


26 


18 


fcrmions 




+ 




I 


26 


18 


fcrmions 






+ 


11+ 






fcrmions 








II- 







TABLE IV. Classification for the doped QDM with resonant 
plaquettes of length 4 and 6. For the square lattice the fam- 
ilies 1+ and I- are equivalent (idem families II, III and IV). 
For the triangular lattice Ji corresponds to the 3 resonant 
plaquettes corresponding to A'^ = 2 in Table II whereas Jg 
corresponds to the second row of A'' = 3 in the same table. 
The two last columns show the references where such models 
have been studied (for Jg = 0) on square (□) and triangular 
(A) lattices. 



of derivation of the QDM Hamiltonian arising from mi- 
croscopic Heisenberg models can be found in Ref. 9 for 
the square lattice with second and third neighbors cou- 
plings and in 10 for the kagome lattice. In these QDM, 
appear a third kind of diagonal or off-diagonal terms 
which has not been considered here, which are dubbed 
assisted terms. An example of such terms is given in 
the last row of table I of 9. They consist of diagonal 
or off-diagonal terms of the kind of the Hv and Hj in 
Hamiltonian (1) but subject to the condition that a third 
dimer is siting in another given neighboring bond. Such 
kind of term can be written as for example : 



(37) 



whose effect is to flip two parallel dimers siting in the 
plaquette i,j,k,l provided that there is one dimer sitting 
in the plaquette m, n. By extending the arguments de- 
veloped above, one can show that under the statistical 
transmutation, this kind of terms transform in the very 
same way as the correspond non-assisted term. For ex- 
ample, the term written above would transform in the 
same way as the term: 

blblihkbu + H.c. (38) 

This is simply due to the fact that assisted terms can 
be written as the product of traditional diagonal or off- 
diagonal operators which we know already how they 
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transform and projectors which are written in terms of 
dimer density operators which are invariant under the 
statistical transmutation. 



V. NUMERICAL INVESTIGATION OF QDM'S 
ON THE TRIANGULAR LATTICE 

A. Summary of phase diagrams in Ref. 18 

We now complement the analytical exact results with 
a numerical study. In our previous work^^, we concen- 
trated on the triangular lattice because it is the best 
laboratory for using our analytic results on the statisti- 
cal transmutation symmetry and for investigating doped 
dimer liquid phases. Here we shall push further these 
studies but we start by a brief summary of the results ob- 
tained in Ref. 18. Considering only flipping terms J = J4 
involving the shortest loops corresponding to A'^ = 2 in 
table (II) a topological (Z2) liquid can be stabilized at 
zero doping^'^^ (the sign of J is irrelevant for a; = 0). At 
finite doping, four non-equivalent families of Hamiltoni- 
ans can be constructed depending on the signs of t and J. 
Note that changing the bare statistics of the holons does 
not introduce a new class of Hamiltonian since this is 
equivalent to change the sign of J as seen in the previous 
Sections. In other words, one can equivalently choose to 
work in the bosonic or fermionic representations. In con- 
trast, the actual statistics of the dressed excitations has 
to be studied numerically. One can use e.g. the method 
developed in Ref. 26 which consists of investigating the 
node content of the wave functions. 

The phase diagrams of the four families of models ob- 
tained in Ref. 18 are reproduced in Fig. 6 for conve- 
nience. At zero doping (where the four models merge 
into the same x — limit) there are two (insulating) 
phases, (i) a six-site cell Valence-Bond-Crystal (VBC) 
for < V/\J\ < 0.7 and a topological Z2 dimer-liquid 
above. At finite doping, family (a) in Fig. 6 is the only 
unfrustrated case and was studied using Green Function 
Monte Carlo (GFMC) methods in Ref. 27. In this model 
bare holons are bosonic and remain representative of the 
physical excitations in the entire region of the phase di- 
agram, as happen also in the Perron-Frobenius square 
lattice version, studied in Ref. 26. The situation is even 
more interesting when J is changed into — J, or equiva- 
lently, bosons are changed into fermions (families (c) and 
(d) in Fig. 6). Finally, changing the signs of both J and t 
of the unfrustrated model (a) leads to the complex case 
(d). Note that the PS regions are further increased at 

< so wc restrict here only to 1^ > 0. 

When J becomes comparable to the holon average ki- 
netic energy (of order xt) holons may be macroscopically 
expelled from the dimer fluctuating background, in order 
to minimize the dimer resonance energy. The question of 
phase separation (PS), i.e. the possibility for the system 
to spontaneously undergo a macroscopic segregation into 
two phases with different hole concentrations, was con- 




FIG. 6. (color online) Qualitative phase diagrams of four 
inequivalent doped QDM's on the triangular lattice versus 
doping (x) and V/\J\ for fixed t/,J = 0.5, from Ref. 18. All 
the models have only fiipping terms corresponding to A'^ = 2 
in table (II). Case (a) corresponds to positive J and t and 
bosonic holes, (b) is obtained from (a) by changing the sign 
of the hoping term t, (c) is obtained from (a) by changing the 
bosons to fermions, and (d) is obtained from (b) by changing 
the bosons to fermions. 



sidered in Ref. 18. In order to perform a Maxwell con- 
struction one can define: s{x) = [e{x) ~ e{0)]/x, where 
e{x) is the energy per site at doping x = n^/N (uh is 
the number of holons in the system and N the number 
of sites). In the case of PS, the energy presents a change 
of curvature at a critical doping x^ corresponding to the 
minimum of s[x) as a function of x. The fact that the 
local curvature of e{x) at a; = is negative then implies 
that the two separated phases will have x = x^ and x = 
(the undopcd insulator) hole concentrations. This study 
revealed that all frustrated (b-d) models show a finite 
PS region (shown in blue in Fig. 6(b-d)) at low doping. 
The extension vs V/\J\ of this region seems to coincide 
with the .T = VBC. In contrast, the non-frustrated (a) 
model does not show phase separation aX V > (as al- 
ready seen in GFMC simulations) but, rather, shows a 
homogeneous region at finite doping where VBC order 
survives. Because of the coexisting non-zero superfluid 
order (U(l) symmetry breaking, see below), this phase 
can be viewed as a "supersolid" (SS). Here super solid- 
ity involves hole pairing in the vicinity of a (insulating) 
VBC phase (and in the absence of PS), as also found in 
the frustrated doped QDM on the square lattice^^ or in 
doped frustrated spin-1/2 quantum magnets. "^^ 

Another important quantity used in Ref. 18 is the sign 
operator defined in Ref. 26 which provides a quantita- 
tive analysis of the nodal structure of the wave function 
and hence gives insights about the statistical nature of 
the holons. i.e. whether they truly behave as bosons or 
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fermions. Such an analysis clearly showed that the GS of 
models (a) and (b) have the same nodal structure as a su- 
perfluid. For models (c) and (d), we dubbed a "complex 
phase" as the statistics of dressed excitations does not 
correspond solely to bosons or fermions. Family (d) also 
shows, in addition, an interesting fermion reconstruction 
at large doping, which we call "fermionic phase" . In this 
work we have extended the study of the nodal structure 
for larger values of the doping than the ones of Fig. 6. 
Our results clearly show that for x > 6 the elementary 
excitation behaves as bosons for the four models. This 
means in particular that in model (d) a dynamical sta- 
tistical transmutation took place in which each Fermion 
bound to a vison in order to form a boson. As we show 
below, this result have important consequences in the na- 
ture of the superfluid phases that we find for those values 
of the doping. 

Finally, an Aharonov-Bohm flux can be inserted in one 
of the hole of the torus, as done in Ref. [29] for the doped 
QDM on the square lattice. A superfluid is characterized 
by well-defined minima in the ground state energy sep- 
arated by a finite barrier in the thermodynamic limit. 
A contrario, a typical signature of (weakly interacting) 
fermions, a flat energy profile is expected even on such 
a small cluster'^". Here, it was reported in Ref. 18 that 
the ground state energy has well-defined minima quan- 
tized at half a fiux quantum for all family of models at 
X ^ 0.25. This might appear as an evidence for condensa- 
tion of chargc-2e particles. However, as already discussed 
in Sec. HI A, the 7r/e-flux periodicity is a generic feature 
of doped QDMs and does not rule out the possibility of 
condensation of deconfined charge-e holons. Next, we 
will characterize more thoroughly the nature of the su- 
perfiuid phases, by introducing a gauge-invariant holon 
Green's function to distinguish the charge-e condensation 
from the usual charge-2e condensation. 



B. New correlations to explore the nature of the 
superfluid phases 

In order to understand the nature of the superfluid 
phases the effective charge of the quasiparticles that con- 
densate have to be determined - either charge-e or charge- 
2e quasiparticles. This is related to the mechanism that 
leads to the spontaneous breaking of the U(l) symmetry 
expected in a superfluid. As a first attempt, one could 
naively use the correlation function (a^aj), but it is not 
compatible with the constraint (8). In other words, it is 
not gauge invariant and thus this correlation function is 
zero. To satisfy the constraint, or equivalently the gauge 
invariance, we need to write correlations in terms of oper- 
ators B. As the hopping of holes can be written in terms 
of operators B we can move one of the holes between two 
distant sites by applying a string of B^ B^s. 

b. Gauge invariant holon Green's function - In the 
subspace where the constraint is satisflcd pairs of holons 
operators a^al acting on sites without holes are equal to 



the identity. Then the holon Green's function we want 
to calculate is given by 

GS^E^^I'^JJa,) (39) 

71 

where 5^"^ = &],„,&ni,«2^«2,n3f'"3,n4-^L-i,n«^"iv.' a 
string operator between the sites i and j following the 
path n. The label (b) indicates that the holes are taken 
as bosons. Similarly, the fermionic version of the Green's 
function is written as 

gS=E(/'4"^/.-) (40) 

n 

where Sj,i = ^],ni^ni,n2^n2,n3^n3.n4---^rijv_i,njv^nN,i- 

Note that, since there arc many ways of moving a hole 
between two sites, the gauge invariance alone docs not 
uniquely determine the definition of the holon Green's 
function. Here we adopt the definition with a sum over 
all possible strings (labeled by n) connecting the two sites 
i and j, with the same coefficient. This definition appears 
most natural to us, as well as in numerical implementa- 
tion. We expect other definitions with some restrictions 
on strings would also work as an order parameter. How- 
ever, the present definition looks advantageous in nu- 
merical calculations, since it can efficiently detect the 
holon condensation with the summation over all possible 
strings. 

Using a complete basis of dimer/hole configurations 
the correlations can be re- written as, 

n a. 13 

= EE(^i")(/5i^)("i// 4"' /^i/^) (42) 

n a. 13 

The holon Green's function can be also written in 
terms of the composite operators only, 

= E^^l'-l,^^^.^ • • • ^2.353,4Sl,,i?l,2) (43) 

G-ij = E(^^v-i,^^^.^ • • • ^LS3.4Sl^^.Bi,2), (44) 

{«} 

where the sum is over all possible paths between sites 
i and j. One can use this representation to show that 
the two correlations are in fact equal, up to an irrelevant 
sign. In other words, we have: 

G'i^] = ±G[f^ , (45) 

where the relative sign depends only on the relative dis- 
tance between i and j. Off-diagonal long range order 
of Gij is a fingerprint of the spontaneous breaking of 
the U(l) gauge symmetry (associated to charge conserva- 
tion). It also implies that the condensing quasi-particles 
have charge e. It is also interesting to observe that at 
large doping, this Green function must necessary have an 
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O bsGrvablGS— ^ 
Phases^ 






(alajaiaj) 




Signs 


SigUF 


Flux periodicity 


PS 


< 


- 


- 


- 


- 


- 


- 


VBC 


> 


LR 


SR 


SR 


- 


- 


2e 


SS 


> 


LR 


LR 


SR 


1 





2e 


2e-SF 


> 


SR 


LR 


SR 


< Sign^ < 1 


< Signj, < 1 


2e 


e-SF 


> 


SR 


LR (weak) 


LR 


1 





2e 


Bose- liquid 


> 


SR 


SR 


SR 


1 





2e 


Fermi-liquid 


> 


SR 


SR 


SR 





1 


2e 


"Complex" phase 


> 


SR 


SR 


SR 


< Signs < 1 


< Sign^ < 1 


2e 



TABLE V. Classification of the possible phases, including various superfluid (SF) phases, that may occur in doped QDM's on 
the triangular lattice. Such phases can be distinguished from the sign of the compressibility re, the long-distance properties 
("SR" means short-range, "LR" means long-range) of various correlations, or the effective charge deduced from periodicity of 
the GS energy versus a magnetic flux inserted through a torus. Sign^ and Sign^ were defined in Ref. 26 to analyze the node 
content of the GS wave function. 



exponential decay (see below). Indeed, the Green func- 
tion for being non zero need at least one path in which 
dimers are present all along it (the string). For a large 
concentration of holes x — >■ 1 it is more and more unlikely 
to find a path with dimers on it so that the Green func- 
tion Gij should roughly decay as (1 — x)^ were L is the 
distance between points i and j. 

Another observable which can detect the charge-e con- 
densation is 

F,,, = {at a]) , (46) 

when a symmetry-breaking groundstate (when it occurs) 
formed by superposition over different dimer-number sec- 
tors is used. Roughly speaking, this corresponds to the 
square of the expectation value of a single hole creation 
operator {a)) in the groundstate, defined in a gauge- 
invariant manner. Such an expression is however less 
convenient to compute numerically in finite-size systems 
and will not be used. 

The scenario of Bosc condensation of polarized spinous 
(the holons in our current formulation) under an applied 
magnetic field advocated in Ref. 27, implicitly implies 
long-range order of Gij. We wish here to substantiate 
such a scenario by an explicit computation of this corre- 
lation function. 

c. Hole pair correlations - If G^j- is short-ranged, 
there is no condensation of charge-e holons. However, 
spontaneous U(l) symmetry breaking and, hence, super- 
fluidity can still occur provided the (hole) pair-pair cor- 
relation, 

P^.,,k,l = {B^.,Bli) (47) 

exhibits long range order. The pair-pair operator is con- 
nected to the square of the single holon effective hopping 



operator in a complicated manner. Namely, 

n n' 

+ («! E'5[^^ E<' 

n n' 

= Bi^jBI I + {loop terms}, (48) 

where the first part of the r.h.s. is obtained from a "clo- 
sure relation" involving all pairs of "retraceable" strings 
n' — n i.e., 

E(4?4?+'5|?45) = fe..4z, (49) 

n 

and the rest corresponds to pair hopping dressed with 
extra loop fluctuations. This suggests that it is physically 
meaningful to write the pair correlations as, 

Pijki = GikGji -f GiiGjk, +Pijki (50) 

where the first two terms can be viewed as the "mean- 
field" contribution and P[ji^i stands for the "connected" 
part in which we remove all the processes involving com- 
positions of single holon hoppings. In particular, both 
sides scale like ((1 — x)"^) when a: — > (a; 1) so that 
it is convenient to normalize Pijki {Gij) by x^(l — x)'^ 
{x{l — x)). Eq. (50) shows that LRO in the holon Green 
function Gik — > Goo characteristic of the charge-e super- 
fluid will induce indirectly LRO in the pair-pair correla- 
tion, P ~ G^. In contrast, the conventional, charge-2e 
superfluid is defined by LRO in the connected part to- 
gether with short-ranged holon Green's function. 

d. dimer-dimer correlations - We finish by recalling 
that the dimer-dimer correlations are expressed in terms 
of the dimer number operators bl jbij as 

A^j,fe,/ = , (51) 

where sites i and j on one hand, and k and I on the other 
hand, are nearest neighbor sites. Long-range order in 
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this correlation function is characteristic of VBC order. 
The wave vector at which the associated structure factor 
diverges defines the VBC wave vector. 

In principle, one can use the new correlations G^.j and 
Pi,j,k,i to refine the previous phase diagrams {Ni^j^^^i was 
used in previous work to determine the VBC and SS re- 
gions). To ease the analysis of the numerical results of 
the doped QDM's, a classification of the various possi- 
ble phases based on simple considerations is provided in 
table V. 

Wc note that, there is no phase where there is a charge- 
e condensation simultaneously with a dimer long-range 
order ("LR" for both {h{jh^jhlihkj) and (aJS'^ja^).) 
This is because existence of the dimer long-range order 
leads to confinement of holons. 



C. Numerical results 




FIG. 7. 16 site-cluster: labeling of the sites (numbered circles) 
and reference bond (purple bond) used respectively in the 
definition of Green's functions and the pair correlations. The 
bonds are labeled according to one of the sites connected to 
them and by a direction, as shown in the example (here site 
10). 

In Fig. 8 are displayed both the holon Green's functions 
(Eqs. (39) and (40)) and the square-root of the pair-pair 
correlations (Eq. (47)) computed by numerical exact di- 
agonalization on a 16-site triangular cluster, varying the 
hole density from x = 0.125 (low hole concentration) to 
X ~ 0.75 (low dimer concentration), from top to bot- 
tom. For convenience, both the Green's functions and 
the square-root of the pair-pair correlations are normal- 
ized by a:;(l — x) to be able to use the same scale for 
all densities. The four hamiltonian classes defined previ- 
ously (see e.g. caption of Fig. 6) are depicted in parallel 
panels (a), (b), (c) and (d) respectively. For all of them, 
we chose the parameters V = 0.3, \ J\ = 1.0 and \t\ = 0.5 
for which, at holon density x = 0.25, the system is either 
in a supcrfluid phase [Fig. 6(a,b)] or in the "complex" 
phase [Fig. 6(c,d)] depending on the QDM class. 



Let us first discuss the data at the lowest hole densi- 
ties X < 0.5. As one can see, only model (a) presents a 
large amplitude of the holon Green's function away from 
the reference site (the largest disk) and its six neighbors. 
While a definite conclusion cannot be drawn from such 
small system, a direct comparison between model (a) and 
the three others reveals a clear change of behavior. In- 
deed, for models (b), (c) and (d), the holon Green's func- 
tion decreases at the largest available distances to signifi- 
cantly smaller values, except maybe for model (c) around 
X ~ 0.5. On the other hand, the pair-pair correlations, 
for which the reference link is in the ui direction, show 
convergence with bond separation to a uniform value for 
model (a) in all relative directions of the two bonds, while 
these correlations are strongly reduced for directions dif- 
fering from that of the reference bond in the other mod- 
els. Hence, our data are clearly compatible with model 
(a) being in the e-superfiuid phase described in Table V, 
unambiguously revealing strong signals simultaneously in 
the holon Green function and the pair-pair correlation. 
Note that we also checked that the dimer-dimer correla- 
tions (not presented here) remain SR in model (a) hence 
reinforcing the previous claim. The behavior of the other 
models at low to moderate doping is less clear, with quite 
smaller amplitudes of Gtj and Pijki at the largest avail- 
able distances. We can however recognize a possible 2e- 
superfiuid phase in model (c) after the phase separation 
zone and up to values of a; ~ 0.2. 

While increasing holon density, from a; = 0.5 to x = 
0.75, we observe that the data for the bosonic and fer- 
monic models become identical, both for t > or t < 0. 
This can be understood by the fact that the (bosonic) 
dimers become then the relevant entities instead of the 
holes. This implies, in particular, that statistical trans- 
mutation or pairing must occur for increasing x for the 
models where fermionic statistics is expected at small 
X (models (c) and (d)). This is indeed confirmed by the 
analysis of the nodal structure of the wave functions that 
we performed for x = 0.625. In this sense the complex 
phase found in model (c) is probably the region in which 
fermions bound to visons and transmute, in order to re- 
semble the bosonic excitations of model (a). Note also 
that charge e superfluidity seems to occur at x = 0.625 
in models (a) and (c) (which seems equivalent for this 
doping). 

However, at larger x corresponding to a dilute gas of 
dimers one expects to eventually recover a charge 2e su- 
perfluid via a continuous (second order) or discontinuous 
(first order) phase transition. This seems to occur al- 
ready for x — 0.75 for models (b) and (d), for which only 
the pair-pair correlations are sizable at the largest dis- 
tances. We have checked that for a lower dimer density 
of 1 — a; ~ 0.1 Gij is short range for all models as reported 
in Fig. 9. In the limit of a very dilute gas of dimers, pair- 
ing between dimers because of the kinetic term J is also 
a possibility. This could result in either phase separation 
or in an homogeneous phase in which both Gij and Pijki 
are short ranged but which is nevertheless supcrfiuid, of 
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FIG. 8. (color online) Holon Green's function (open circles) and square root of the absolute value of the pair-pair correlations 
(filled symbols) for parameters V = 0.3, \ J\ — 1.0 and |i| = 0.5, at various densities ranging from x — 0.125 to s = 0.75 (2 to 
12 holons on 16 sites) and for the 4 classes (a-d) of models defined in the caption of Fig. 6. The pair-pair correlation Pi,j,k,i is 
defined by a reference bond orientation fj — fi = Ua and the orientation of the final bond ti — Tk = up. In our case, we chose 
Ua = Ml and we consider three cases for up: ui (filled squares), U2 (filled circles) and ug (filled diamonds) - see Fig. 7. 
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FIG. 10. Phase diagrams of the four models (a), (b), (c) and 
(d) at 1^/|J| = 0.3 and t/\J\ = 0.5 derived from Fig. 6 and 
Fig.8. All have both the e-superfluid (e-SF) and 2e-superfluid 
(2e-SF) phases at different doping depending on the model. 



phase followed by a charge 4e superfluid phase in models 
(b) and (d) due to dimer pairing. 



FIG. 9. Holon Green's function (open circles) and square 
root of the absolute value of the pair-pair correlations (filled 
symbols) for parameters V = 0.3, \J\ = 1.0 and \t\ = 0.5, 
at a low dimer density 1 — x ~ 0.1 for the 4 classes (a-d) of 
models defined in the caption of Fig. 6. Left; 2 dimers on a 
36-site cluster (x ~ 0.89). Right: 3 dimers on a 64-site cluster 
{x ~ 0.91). 



coherent dimer pairs of charge 4e. We have checked that 
there is not phase separation in none of the four mod- 
els at those large values of x. We have also looked at 
the energy difference between two and one dimers and 
found that pairing is indeed favored in models (b) and 
(d) which may explain the drop of the Pijki in Figure 9 
for X > 0.9. 

Based on Figs. 6, 8 and 9, we have extracted the 
qualitative phase diagrams for the four models at fixed 
V/\J\ — 0.3 and t/\J\ = 0.5 as a function of doping 
X. They are depicted in Fig. 10. Charge e superfluidity 
seems to occur in all models, with the largest occurrence 
in model (a). For intermediate doping, models (b) and 
(d) seem to present short range correlations for both one 
and two particle Green functions. This behavior sug- 
gest an uncondensed phase which in the case of model 
(d) would correspond to a Fermi-liquid state. Since el- 
ementary excitation in model (b) are bosonic the pres- 
ence of an uncondensed phase points toward an exotic 
Bose-liquid state, although this statement should require 
a more detailed study (using clusters of a much bigger 
size) which is beyond the scope of the present article. 
For X big, as expected, all models exhibit a 2e-supcrfluid 



VI. CONNECTION TO BOSE-HUBBARD 
MODELS 

We finish this work by discussing the connection to 
Bosc-Hubbard like models which do not contain a priori 
the ice-rule constraint. However, the physics of the doped 
QDM can emerge naturally when some form of large re- 
pulsion between the itinerant bosons is considered, hence 
providing emergence of fractionalized excitations.'^^ 

In Ref. 32 it was introduced a simple model of hard- 
core bosons hopping (t) on a kagomc lattice with a boson 
repulsion favoring the smallest number of bosons in 
each hexagon. 

(ij> O 

where d\ creates a boson on site i and Uq = X]i=i d\di 
is the number of bosons in a hexagonal plaquette. When 
the boson density is p = 1/2 a large Vq/I stabilizes an in- 
sulating phase whose quantum dynamics is described by 
a generalized QDM on the triangular lattice with exactly 
3 dimers per site. The insulating phase is a Z2 topological 
liquid. When t < the model is not frustrated and can 
be studied with QMC : the superfluid-insulator transi- 
tion was argued to be a novel non-conventional fractional 
critical point.^^ 

To make the connection with some of our doped 
QDM's we shall assume here the microscopic d bosons 
have charge — 2e and their density is set to p = ^(1— x/2), 
x << 1. In that shown in Fig. 11, for Vo/t — !■ 00 

the lowest-energy configuration space {E ~ A^V^/3) is 
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FIG. 11. (color online) Hard-core bosons on the kagome lat- 
tice. Covering the lattice with one, and only one, boson per 
hexagonal plaquette (ice-rule constraint), we can identify dis- 
tribution of bosons with a dimer hard-core covering on the 
triangular lattice. Removing a boson of charge 2e creates two 
defect hexagons (shaded). Each defect (hole) has charge e and 
can move on the lattice. Coherent hoping of two hard-core 
bosons corresponds to a flipping process in the dimer model. 



given by all hardcore dimer- coverings on an effective tri- 
angular lattice, where each d-boson has been replaced 
by a dimer connecting two sites of the triangular lat- 
tice. Such configurations respect a local ice-rule con- 
straint with one, and only one, boson per hexagon. When 
X = 0, moving a single d-boson violates this ice-rule so 
one has to move at least two simultaneously. This pro- 
cess of amplitude J ~ t'^/^j corresponds exactly to a 
dimer flip on a lozenge, identical to the one of the QDM. 
Strictly speaking this mapping onto the QDM does not 
involve any dimer-dimer repulsion V. However, one can 
add a small third-nearest neighbor density-density repul- 
sion Vdd << ^ between the d-bosons located on differ- 
ent hexagon. In the mapping for large T^/i, this inter- 
action translates directly into the dimer-dimer repulsion 
V — Vdd- Thus, by tuning Vdd in the Bose-Hubbard 
model, the topological Z2 (insulating) liquid can be sta- 
bilized. 

When X ^ an empty site on the original kagome lat- 
tice, corresponds to two "defect" hexagons carrying an 
overall charge 2e w.r.t. the insulating GS. It is easy to see 
(Fig. 11) that each "defect" hexagon can move indepen- 
dently on the effective triangular lattice by simple pro- 
cesses that involve a single d-boson hopping. Therefore, 
each "defect" can be considered as an effective charge-e 
hole on the triangular lattice. The amplitude t of the 
hole hopping is the same as the one of the microscopic 
d-boson Hubbard model. When a d-boson of charge — 2e 
hops by a lattice spacing a, the effective hole of charge e 
hops by a distance 2a so that the charge center of mass is 
conserved. Note also that the hole density in the effective 
QDM is X. 



FIG. 12. (color online) Schematic and speculative phase di- 
agram of interacting HCB on the kagome lattice for Vdd = 

(a) , with fine tuning of the third- nearest neighbor repulsion 

(b) - see text. In (a), the exact nature of the transitions be- 
tween the 2e-superconductor and the VBC insulator at a; = 
needs to be further investigated. In (b), the dot at a; = 
might correspond to the XY* fractional critical point be- 
tween the superfluid and the topological Z2 insulator. The e- 
superfluid is the same phase as in the doped QDM of Fig. 6(a). 



The mapping to the doped QDM on the triangular 
lattice is therefore complete. However, it is important to 
notice that J = f^/VQ > (for real t) so that only QDM 
(a) and (b) can be realized with HCB depending on the 
sign of t. Introducing imaginary hoping t ~ ir on the 
kagome lattice equivalent to put U(l) fluxes through the 
triangles lead to the QDM models (c) and (d) in the pres- 
ence of a magnetic field. In practice only the case of a real 
t < hopping on the original kagome lattice can be han- 
dled with QMC. Assuming the phase of the correspond- 
ing doped QDM ((a) model) is a fractionalized charge-e 
superfluid, we therefore predict a non-conventional phase 
transition between the (ordinary) 2e-superfluid of the 
weakly interaction d-bosons and an exotic e-superfluid 
at large V^/i, as schematically shown in Fig. 12. This 
is possible if the third-nearest neighbor repulsion is care- 
fully tuned ~ Vdd ^ J = /Vq. If Vdd = 0, one gets 
a transition to a plaquette VBC phase at a; = 0, which 
might involve intermediate phases. Indeed, close simi- 
larities are expected with the melting of the (bosonic) 
plaquette VBC on the checkerboard lattice, revealing 
an intermediate commensurate supersolid.'^'^ Similarly to 
the effective triangular QDM at y = 0, doping of the 
VBC insulator should immediately result into a super- 
solid phase which would melt into a charge-e superfluid 
above some critical doping. Lastly, at even larger doping 
(corresponding to a dilute gas of dimers) a second phase 
transition to a charge-2e superfluid is expected. 
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VII. CONCLUSIONS 

In this paper we have estabhshed a rigorous and 
general equivalence between QDM Hamiltonians with 
bosonic holes and a corresponding QDM Hamiltonian 
with fcrmionic holes. Although this correspondence was 
already noticed on the basis of numerical simulations in 
Ref. 26 and established analytically in Ref. 18 for Hamil- 
tonians with the simplest flipping term, the correspon- 
dence has now been generalized to more complicated 
cases. More importantly, we provide a general recipe 
to very quickly - and without any computation - estab- 
lish which are the two equivalent Hamiltonians under this 
statistical transmutation. 

We also note that, when working with finite size sys- 
tems, while the composite particle representation is valid 
for any kind of boundary conditions, this is not the case 
for the method that use the Jordan- Wigner transforma- 
tion. Indeed, the issue of boundary conditions in the 
two-dimensional version of the Jordan- Wigner transfor- 
mation has been very little discussed in the literature. 
The point is that it does not seem possible to impose 
periodic boundary conditions in a consistent way when 
using the two-dimensional version of the Jordan- Wigner 
transformation, even if the total numbers of particles 
is kept fixed. This can be contrasted with the one- 
dimensional version of the transformation where periodic 
boundary conditions can be consistently imposed pro- 
vided one keeps the number of particles fixed. In this 
sense the analytical results obtaincd^^ with the help of 
the Jordan- Wigner transformation are only valid to in- 
finite systems or finite system with open boundary con- 
ditions while the composite particle representation used 
here can be consistently applied for any kind of bound- 
ary conditions. We have then provided many examples 
of equivalent Hamiltonians for the more generic cases of 
the square, triangular, hexagonal and kagome lattices. 

In order to characterize the exotic superfluid phase 
with condensation of fractionalized charge-e holons, we 
have introduced the gauge-invariant holon Green's func- 
tion. We have then considered four inequivalent cases 
of QDM Hamiltonians on the triangular lattices and 
have numerically studied various correlation functions in- 
cluding the above mentioned holon Green's function, by 
Lanczos exact diagonalization of finite-size clusters. We 
obtained rather strong and direct evidence for the exis- 
tence of the exotic superfluid phase due to condensation 
of holons carrying charge e. In fact, our numerical results 
suggest that all the four models we have studied exhibit 
the charge-e superfluid phase. More conventional chargc- 
2e and 4e superfluid phases are also present. While the 
existence of a charge-e superfluid phase may be natu- 
rally understood when the bare particle are bosons, it is 
much more puzzling in the case where they are fermions, 
as it corresponds to a superconductor without Cooper 
pairing. Although one may expect that the condensa- 
tion of charge e holons requires the holes to be bosons, 
our rigorous mapping shows that fcrmionic statistics can 



be always assigned to holes in the microscopic Hamil- 
tonian. This implies a dynamical statistical transmuta- 
tion of holons in the QDM where holes are represented 
as fermions. Those kind of dynamical statistical trans- 
mutations can be monitored by studying the nodes of 
the wave function and one spectacular example can be 
found in model (d) at intermediate values of the doping 
where the system seems to switch from a Fermi liquid to 
a (bosonic) e-superfluid phase. 

We have then provided a concrete microscopic Bose- 
Hubbard Hamiltonian which in the strong interaction 
limit behaves as a QDM on the triangular lattice, as the 
one analyzed numerically. It allow us to have a better 
control on the doping (by simply varying the number of 
bosons) and to better visualize the superfluid phases e- 
superfluid and 2e-superfiuid. It is important to stress 
that a QDM may arise as effective low energy models of 
quite different microscopic Hamiltonians. As such, the 
physical consequences of various phases in the QDM can 
depend on the mapping. For example, if the QDM arises 
from a microscopic electronic Hubbard model, the holes 
are real electrons vacancies (models (c) and (d)), and 
then they are charged. In this case the different super- 
fiuid phases are superconducting phases. There is how- 
ever another way in which one could introduce doping. 
Imagine for example a system in which there is no real 
electron vacancies but some magnetic field is applied to 
the system. The effect of the magnetic field may have as 
effect to break some of the singlet that are represented 
by the dimers leaving two polarized spin 1/2, which now 
play the role of the holes^^ (models (a) and (b)). In this 
case the holes are neutral but carry spin, so that the su- 
pcrfiuid phase now corresponds to a superconductor of 
magnetic current. Our results provide a validation of the 
previous claim^^ of an exotic superfluid of condensed de- 
conflned and polarized spinous (equivalent to our charge 
e-superfiuid). In addition, we predict here the existence 
of another phase of deconfined spinous, the Bose liquid, 
corresponding to an exotic spin liquid carrying uncon- 
densed (polarized) spinous. Interestingly, such exotic 
phases could indeed be realized in simple frustrated mag- 
netic systems, as for example the kagome anisotropic spin 
1/2 model close to the magnetization plateau at 1/3 of 
the saturation value. We hope that the results of this 
paper will establish new motivations to investigate, with 
a new light, different microscopic models which may give 
rise to the doped QDM as an effective low energy model. 
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Appendix A: Connection with the two-dimensional 
Jordan- Wigner transformation 

In this appendix we elaborate on an alternative proof 
of the statistical transmutation in the QDM. In this ap- 
proach we use a two dimensional version of the Jordan- 
Wigner transformation. The fundamental ingredients of 
this transmutations were presented for the square and 
triangular lattices in a previous work"'^^. It is important 
to stress that this procedure is totally generic and can 
be implemented in any two-dimensional lattice with open 
boundary conditions^^'^^. As the main technical steps for 

I 



Hj=Je 



Hv=Ve 



.A' 



(b) 



Ac) 



the square and triangular lattice were already presented 
in 18, here we only show the details for the kagome lat- 
tice. As we have discussed above, for this lattice there is 
a big freedom in choosing the ordering prescription. This 
fact given rise to an extra freedom in taking the sign of 
the flipping constants when we change the statistics of 
holes. Wc will see in this section that this freedom is ma- 
terialized whithin the Jordan- Wigner approach by using 
gauge transformations on dimers and holes. Let us start 
with a quantum hard-core dimer model in the presence 
of holons on the Kagome lattice given by the following 
Hamiltonian; 



H = Hj + Hv + Ht 



(Al) 



The terms Hy and Hj corresponding to the diagonal and 
off-diagonal terms of the pure dimer model are taken up 
to resonance plaquettes of length 1223.24^ 



^A=, 




(A2) 
(A3) 



The sum runs over all the hexagons of the lattice and all 
the possible orientations of the plaquettes are implicit. 
The kinetic and diagonal terms are given by 



H.c. 



(A4) 



•rT-\/TT 



and similar expressions for the other terms. For conve- 
nience, the labels a in the amplitudes Jq, and Va cor- 
respond to the length of the associated resonance loops 
and when it corresponds we add the label (a), (6) or (c) 
corresponding to the 3 non-equivalent plaquettes for the 
cases of length 8 and 10. The Hamiltonian corresponding 
to the hoping of holons is given by 



Ht = H 



it) 



it) 



(A5) 



where 



H 



it) 

A 



A 



{|a)(o\| + |/o)(a| + |o\>|/ 



H.c.}. 



(A6) 



The hoping of holons can be written in a general way, 
independently of the lattice, as a sum of three-site Hamil- 
tonians 



(ijk) 



(A7) 



with 



\ijk) 



(A8) 



Where, we have been projected the Hamiltonian on the 
subspace where the constraint 



(A9) 



is satisfied by mean the projectors V. Where the sum 
runs over NN of site i. Starting from h^^^^^ we trans- 
form the boson operators ai using 



(AlO) 



with 



'^' = E/l/^^'^g(^j-^') (All) 

togheter with the following transformation for the dimer 
operators 



1-] l;3 



I. . — pii<Pi+<t>j) U. . 



we obtain 



(A12) 
(A13) 

(A14) 
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where the hoping amphtude is given by 



(A15) 



Equations (A14) and (A15) are, in fact, independent 
of the lattice details. The information concerning the 
lattice geometry is contained in the arguments on the 
exponential of Eq. (A15). This equation can be written 
in a compact form as 



t = -te 



(A16) 



where tpt = arg(Tj — r^) — arg(Tj — Tk) can be represented 
graphically as 



K K 

A, -A, 



(A17) 



K 

and represents the aTg{Tj ~ n). 

i J 

Now, we study the kinetic Hamiltonian corresponding 
to dimers. Let us start with the smallest resonance loop 
compatible with NN dimers on the Kagome lattice, the 
plaquete of length 6. In this plaquette, the resonance of 
the two possible dimerizations is given by 



H 



Je 



H.c. 



(A18) 



In order to transform the dimers to the new representa- 
tion using the flux generated by the statistical transfor- 
mation of the holons we write the Hamiltonian in terms 
of dimcr operators bi,j as 



(A19) 



with 



{i,j,k,l.m,n) 



Je ^lA ihl,^J},,kbL„^bn^; + H.c.(A20) 



Using the transformation (A12) is straightforward to 
write the Hamiltonian as 



where 




= ( arg(fm - r„) + arg(fi - fj) + arg(ffc - n) 



((arg(r„ 



Ti) + arg(fj - Tk) + arg(fi - t, 



The two graphs correspond to the initial and final dimer- 
ization on the plaquette. In each graph we replace the 
dimers by arrows drawn in a clockwise direction and each 
graph represent the sum on the arguments of the arrows. 

For the resonant plaquettes of length 8 we have 3 topo- 
logically distinct configurations. Let us study now the 



term corresponding to the resonance plaquette 



After write it in terms of dimer operators and transform 
following (A12) we obtain 



J. 



(&) 



with 



V'; 



(fe) 




(A24) 



(A25) 



where 




= (^arg(f,„ - f„) + arg(ffe - fj) + arg(fi - fj) 
+ arg(fo - Tp) ) 




- To) + arg(fi - fvi) -I- arg(fj - ft) 



+ arg(fp - Ti 



{i,j,k,l.rii,n) 



T 7i 7i 7i 7 7 7 tx / a r.-. \ It is casv to chcck that, under transformation (Al2), 

^6 blAA,nhkbLr.A., + H.C.(A21) ,(7) L ' 



where 
and 

-06 = 



(A22) 



(^aTg{Trn - T„) + arg(fj -~fj) + arg(ffc - ti)^ 

- ((arg(T„ - Ti) + arg(fj - ffc) + arg(fi - f^)^ . 

It is convenient to use a graphical representation for the 
phase ipQ 



V'6 




(A23) 



the amplitudes jj^^ corresponding to the rest of the res- 
onance plaquettes also transform as 

J_('^) = -4'')e''>^"', (A26) 

where ip^^ is the phase obtained from the difference be- 
tween the two possible dimerizations in a given resonance 
plaquette of the sum of the arguments corresponding to 
dimers oriented clockwise (or anticlockwise). 

This graphical rule can be used to study higher or- 
der therms in the kinetic Hamiltonian. This allows to 
determine the Hamiltonian after the J-W transformation 
on the holons. In the Kagome lattice, up to resonance 
plaquettes of length 12 we obtain that the Hamiltonian 
Hj corresponding to kinetic energy of the dimers can be 
written as in (A2), but dimers arc now created by the 
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TABLE VL Values of Jc,/ Ja corresponding to the lowest or- 
ders of the resonant plaquettes. 



operators j and the amplitudes Ja must be replaced 

by Ja where as the values of Va remains unchanged. In 
table VI we show the values Ja/ Ja for the resonance pla- 
quettes up to length 12. 

After the transmutation, the Hamiltonian correspond- 
ing to the hoping of holes is given by 



Ht = H 



(t) 



H 



it) 



(A27) 



where in dimers are also created by operators 5| , 



and the holes are fermions created by the operators 



Hi'' = i 



E{|a>(o\| + |/ 



H.c.}, 



A) / 

(A28) 



where i = te^'a'^. In order to obtain the original value 
of the hoping constant for the holons {t — )- t), we can 
perform a gauge transformation on the dimers. This is 
possible, for example, by mean the following gauge trans- 
formation 



independently if the dimers are on up or down triangles. 
We have for the hole hoping term 



H.c.} 



A) ! 

(A29) 



while for the resonance terms the gauge transformation 
does not change the couplings J and Y . Then we obtain 



U, = J, 




Finally we can use another gauge transformation to get 
a more simple Hamiltonian. We change the dimers on up 
triangles as 



A)^^'''^ a; 

while we change the corresponding to down triangles as 

v)^e-|v; 



Taking the values d\ — — H^i — '^2 = 0, 6*3 = 7r/2 and 
(y93 = —'it/2. The amplitudes J a change as 



Hj = -Je 




Then finally we obtain a Hamiltonian with a global 
change of sign in the amplitudes of the kinetic term. Af- 
ter the transformation, the terms corresponding to the 
hoping of holons becomes 



/2 



+ \A 

+ H.c. 



^>|/o>} 



(A30) 



This hoping Hamiltonian can be transformed by mean a 
simple gauge transformation on the holes, fj — e^'^"^\fj 
with Q = ( f 1 ) I to recover the original form of the 
Hamiltonian. With this transformation we have that 

and we recover the original form (A29) for the Hamilto- 
nian corresponding to the holes. Then, finally we have 
that changing the statistics of holons together with the 
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sign of all the kinetic amplitudes we obtain a completely 
equivalent Hamiltonian. 

By using the J-W transformation we can recover the 
equivalences presented in previous section, but the pro- 
cedure is more laborious. We have pointed that starting 
from different prescriptions for the bonds we can prove 
different equivalences. A change i the lattice prescription 
used in the composite operator representation is equiva- 
lent to a gauge transformation on on the dimers. 

One last very important point concerns the issue of 
boundary conditions within the J-W approach. The com- 
posite operator approach which we have extensively used 
in the paper is valid independently of the boundary con- 
ditions used for the system. In contrast, the J-W ap- 
proach is valid only for infinite systems or finite systems 
with open boundary conditions. Let us come back to 
equation (AlO) assuming that the bosonic a operators 
are well defined in a system in periodic boundary condi- 
tions. This means that for example the operators and 



ai are forced to be the same if sites i and i correspond 
to the same point in the system because of the periodic 
boundary conditions. However, because of the very na- 
ture of the non-local J-W transformation, the relation 
between the operators /,; and fi must contain a twist. In 
the one-dimensional version of the J-W transformation, it 
is easy to see that this twist is a sign which depends only 
on the total number of particles in the system. Then, 
in one dimension, if one restricts to the subspace of a 
fixed number of particles it is possible to implement con- 
sistent way periodic boundary conditions. Here in two 
dimensions, the twist one should force for the fermionic 
operators depends not only on the number of particles 
but also on their relative positions with respect to the 
points i and i. That means that, even restricting to a 
fixed number of particles, starting from periodic bound- 
ary conditions for the bosons operators, it is not possible 
to implement boundary conditions for the fermionic op- 
erators which are consistent with all the possible particle 
configurations. 
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